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ABSTRACT 

We present new constraints on the merger history of the most X-ray luminous cluster 
of galaxies, RXJ1347. 5-1145, based on its unique multiwavelength morphology. Our X-ray 
analysis confirms that the core gas is undergoing "sloshing" resulting from a prior, large-scale, 
gravitational perturbation. In combination with multiwavelength observations, the sloshing gas 
points to the primary and secondary clusters having had at least two prior strong gravitational 
interactions. The evidence supports a model in which the secondary subcluster with mass 
M=4.8±2.4x 10^^ Mq has previously (>0.6 Gyr ago) passed by the primary cluster, and has 
now returned for a subsequent crossing where the subcluster's gas has been completely stripped 
from its dark matter halo. RXJ1347 is a prime example of how core gas sloshing may be used 
to constrain the merger histories of galaxy clusters through multiwavelength analyses. 
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1. Introduction 

Clusters of galaxies represent the end state 
in the hierarchical growth of the largest gravita- 
tionally bound structures in the universe. Out- 
side their cores, the (near) universal temperature 
and mass profiles of relaxed clusters make them 
ideal tools for constraining cosmological parame- 
ters (e.g. wq and as) whic h govern the growth o f 



structure in the universe ( Vikhlinin et al. 20091 ) 



We see this growth of structure in l arge scale sim- 



ulations (e.g. Millennium I and I I; ISpringel et al 



20051 : iBovlan-Kolchin et al.ll2009l ). however, there 



are relatively few corresponding observational con- 
straints in place. This is particularly true for 
matching the predicted rates at which massive 
clusters assemble with direct observational con- 
straints on those rates. For individual clusters, 
one should use empirically derived initial condi- 
tions in order to constrain their merger histo- 
ries and attempt to reproduce the observations in 
han d, as was done recen tly for the Virgo Clus- 
ter (jRoediger et al. I I2OIII) . " One such cluster for 



which there exists a wealth of observational data 
is the most luminous X-ray cluster, RXJ1347.5- 
1145 (z=0.451). 

High resolution hydrodynamic simulations pre- 
dict the effects that multiple gravitational inter- 
actions will have on galaxy clusters. In particu- 
lar, the resulting morphology of the low-entrop; 



iar, tne resulting morpnology or tne low-entropy 
cluster core gas ( Ascasibar fc MarkevitchI 2006t 



Poole et al.ll2008l : IZuHone et al.ll2010l : lRoediger etal 



20111 ). In this paper, we use the multiwavelength 
observations, along with simulations, to propose 
a merger history for RXJ1347 which explains the 
current morphology of the cluster in optical. X-ray, 
and millimeter Sunyaev-Zel'dovich (SZ) observa- 
tions. In the following section, we describe the 
phenomenon of core gas sloshing in galaxy cluster 
cores. In Sjll we present a brief discussion of the 
extensive multiwavelength coverage and analysis 
of the cluster. In f|3]we discuss our analysis of the 
Chandra X-ray observations and the identification 
of the sloshing gas in the cluster core. Finally, 
in Sj4] we present new constraints on the merger 
history of the cluster based on the accumulated 
multiwavelength observations. 

1.1. Core Gas Sloshing 

The model for core gas sloshin g has been ex- 

plored numerically in several works (ITittlev fc Henriksen 



20051:1 Ascasibar fc Markevitchl2006l:IZuHone etal 



2OIOI : iRoe digcr et al. 2011j), and id entified obser 



vatiqnally in many galaxy groups (IRandall et al 



2009; 



2001 



2007 



Machacek et al. 



Churazov et al 



2010') an d clusters (iMarkevitch et al 



2003; Markevitch fc Vikhlinin 



Owers et al.ll2009HJohnson et al.ll2o'To[ r 



The sloshing phenomenon begins when two 
group or cluster-sized halos cross nearby one an- 
other (<1 Mpc depending on their masses). As 
seen from the reference frame of the more massive 
object (the primary cluster), the less massive ob- 
ject (the subcluster) falls in on either an elliptical 



or hyperbolic orbit (the exact orbit is irrelevant 
for the subcluster's first crossing). As the subclus- 
ter falls in, the primary cluster (its dark matter 
halo and cluster gas) experiences a gravitational 
force towards the incoming object. The subcluster 
reaches pericenter in its orbit, at which point the 
gravitational force is at its peak, and after which 
point the primary cluster reaches its maximum 
displacement from its initial position. 

The wake of the passing subcluster directs the 
ambient velocity field towards the primary clus- 
ter, while at the same time the primary clus- 
ter continues to move towards the subcluster. 
This leads to ram pressure on the primary clus- 
ter gas, causing it to compress along the lead- 
ing edge of its trajectory. As the subcluster 
recedes, the velocity field around the primary 
cluster changes rapidly, removing the ram pres- 
sure, and resulting in a ram pressure "slingshot" 
( Hallman fc Markevitch|[2004 ) , where the primary 
cluster gas moves quickly ahead of the potential 
minimum, causing a separation between the two. 
It is this initial separation that leads to the sub- 
sequent production of large-scale asymmetries in 
the gas distribution. 

After its initial separation, the primary clus- 
ter gas density peak (its "cool core") falls back 
towards the potential minimum, which has itself 
been displaced from its initial position. The den- 
sity peak moves against the fiow of its own trail- 
ing gas, resulting in ram pressure along its lead- 
ing edge. This ram pressure causes a significant 
(> 2x) compression of the gas along the edge, 
revealing itself as a crescent-shaped edge in the 
X-ray surface brightness distribution of the core 
gas. Upon reaching the apocenter of its oscilla- 
tion about the potential minimum, the cool core 
reverses direction and again falls back towards 
the potential minimum against the flow of its 
own trailing gas, resulting in the formation of an- 
other ram pressure fro nt. These fronts have been 
jiubbed "cold fronts" (| Vikhlinin et al.ll2Q0ll) . as 
opposed to shock fronts, because they delineate 
the edge of the cool core's oscillatory motion, with 
the coolest gas interior to the front and the warmer 
gas outside. This oscillation and front production 
continues for long periods (>1 Gyr) in the core 
but is progressively damped with each crossing of 
the density peak through the potential minimum. 

The above describes the essential components 
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of the radial motion of the density peak about 
the potential minimum. However, because angu- 
lar momentum is transferred from the subcluster 
in this interaction, the density peak undergoes az- 
imuthal, as well as radial, oscillation about the 
potential minimum. In all but a head-on merger 
(&=0), the displaced gas core does not pass directly 
through the potential minimum, but rather un- 
dergoes a damped oscillation around the potential 
minimum. Whether due to a head-on or off-axis 
collision, the resulting oscillatory motion of the 
lowest entropy core gas within the gravitational 
poten tial is called "sloshing" (jMarkevitch et al 

2nnih . 



Despite the dependence of the long-term evo- 
lution of these cold fronts on merger parameters 
(cluster mass ratios, impact parameters, etc.) and 
the plasma properties (e.g., viscosities, magnetic 
fields, etc.), there are several qualities of these 
sloshing cold fronts that appear universal and are 
most relevant to RXJ1347. First is that, based on 
the simulations, there appears to be a minimum 
time after subcluster pericenter for the first cold 
fronts to form. This is because the dense gas core 
in the primary cluster must first slingshot past the 
cluster's potential minimum and will only then 
form the first cold front as it falls back towards 
the potential minimum. Though the ability to de- 
tect the first sloshing cold front formed is limited 
by the simulation spatial resolution, over an order 
of magnitude in both mass ratio and impact pa- 
rameters, simulations show the first cold front does 
not form until ftiO. S Gyr after the closest approach 



p proacn 
hi 120061: 



of the subcluster (lAscasibar fc Markevitcf 

ZuHone et cH l20ld iRoediger et al.l l201l[ ). Sec- 
ond, the location of the cold fronts reveals when, 
and from which direction, the sloshing disturbance 
originated. The radial distance (projected or not) 
from the X-ray peak where the cold fronts are 
found always increases with the time since cold 
front formation. That is to say that the sloshing 
cold fronts which formed first are always further 
from the oscillating core than fronts which form 
at a later time. Also, simulations show that the 
first cold front always forms on the opposite side 
of the cluster as the location of the initial distur- 
bance. This is because the initial displacement of 
the primary cluster's gas core occurs on the same 
side as the pericenter of the disturbing object with 
the first cold front then forming as the core falls 



back from this initial displacement. And finally, 
simulations show that the inclination of the clus- 
ters' merger axis with respect to our line of sight 
dictates whether the cold fronts will appear to us 
as a spiral inflow (for a merger in the plane of 
the sky), concentric crescent shapes (for a merger 
along the line of sight), or some combination of 
the two resulting from an inclined spiral inflow. 
We discuss the implications of these for the his- 
tory of RXJ1347 in 



2. Previous Observations of RXJ1347 

As the most X-ray luminous cluster known 
(with Ly=6.0 X lO'"^ erg s^^ from 2 to 10 keV; 
Voges et al. I lf999l) . RXJI347 has been extensively 
studied in the radio, millimeter, s ubmillimeter 
(|Gitti et al.ll2007al:lMason et al.ll2010l). in the opti 



cal via both weak and strong lensing ([Bradac et al 
20081 : iMiranda et al] 120081 ) . as well as through 



spect r oscopy of memb er galaxies ([Cohen fc Kneib 
l20 02t ILu et a l.' '20I0 l), and of course in X-rays 



(|Allen et al.l i2002; .Gitti et al.l l2007bl) . Indeed 



one could argue that of those clusters beyond 
=0.3, only the "Bull et" cluster (IE 0657-56; 



Markevitch et al.l (|2004f ) ) has received more atten- 



tion. To establish the context and relevance of our 
interpretation, we review here some of the previ- 
ous observations of RXJI347. We focus primarily 
on those contemporary works most relevant to our 
discussion. 

2.1. Previous X-Ray Observations of RXJ1347 

RXJ1347 (see Fi g. □ was first observe d in X- 
rays with ROSAT (ISchindler et"aDll995h . With 
the arcminute angular resolution of the ROSAT 
PSPC, RXJI347 appeared to be a relaxed sys- 
tem with a central cooling peak. Since ROSAT, 
the cluster has been observed b y nearly every X - 
ray satellite, in cluding Suza ku (lota et al.ll2008[ ). 
XM M-Newton (iGitti et al.l l2007bl) and Chan- 
dra (| Allen et all 120021: iBradac et al.ll2008l ). Con- 
clusions regarding both the overall gas density, 
temperature, and metallicity are in agreement. 
RXJ1347 exhibits a hot {kT >I0 keV) intra- 
cluster medium (ICM), with an extension to the 
southeast. Prior to the Chandra observations 
(jAllen et al.l |2002| ). there was little indication 



of this southeast structure (RXJI347-SE here- 
after) and, excluding this quadrant, the cluster 
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Fig. 1. — Combined 73 ks Chandra image of RXJ1347.5-1145 from OBSIDs 506, 507, and 3592. The image is shown in the 
0.5-2.5 keV band and has been exposure corrected and background subtracted. The surface brightness edges can be seen on 
either side of the emission peak. The two small circles denote the locations of the primary and secondary cD galaxies (cDl and 
cD2 respectively) and the small blue circle denotes the location of the X-ray peak. In all images, north is up and east is to the 
left. A color version of this figure is available in the online journal. 



gas distribution appears s ymmetric ou t to la rge 
scales (>1 Mpc), leading lAllen et all (|2002[) to 
conclude that "the cluster appears relatively re- 
laxed." High-resolution imaging and spectral ob- 
servations, with both XMM-Newton and Chan- 
dra, reveal a more complex picture with several 
features not typically seen in relaxed clusters in 
hydrostatic equilibrium, including a non-isotropic 
temperature distribution and two extremely hot 
regions (fcT > 20 keV) within 100 h'^ kpc of the 
X-ray brightness peak (see Fig. [2fc)). 

2.2. Optical Observations of RXJ1347 
Cohen fc KneibI (|2002[ ) and lBradac et al.l (|2008[) 



show that the cluster hosts two cD galaxies, one 
0.5" (3.3 h-^ kpc) NE of the X-ray peak (labeled 
cDl in Fig. ID) and the other, ~ 18" (118 h'^ kpc) 
to the east, which has no detected X-ray counter- 
part (labeled cD2 in Fig.[T]). The recessional veloc- 
ities for these two galaxies differ by only ^ 200 km 



s-i (|Cohen fc KneibI 120021: [Lu et al.ll2010l ) which 
is much less than even the lowest estimate for the 
velo city dispersion of the c luster (cru=910±130 km 



Cohen fc KneibI |2002[ ). 



A dedicated study was performed using 400 
redshifts in the region surrounding RXJ1347 



(|Lu et al.l I2OIO ) which, alo n g wit h the previous 



efforts by 
number of cluster 
shifts to 



Cohen fc KneibI (|2002l ). brought the 
members confirmed via red- 



140 galaxies. iLu et al. I (l2010h note sev- 
eral new structures in velocity space, located co- 
linearly in projection along a roughly NE /SW line 
throu gh the cluster center (see Fig. 1 in lLu et al.l 
2OIO1) . Among these galaxy overdensities is a sig- 



nificant {7a) peak in the redshift distribution cor- 
responding to a cluster of galaxies approximately 
~7 h"^ Mpc (projected) SW of the main cluster 
(referred to as RXJ1347-SW hereafter). Ahhough 
this structure appears to be much more distant 
than RXJ1347 (its mean redshift is ~4000 km s"^ 
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larger), there is evidence that it is physically as- 
sociated with RXJ1347, either through a previous 
interaction, which is unlikely given the separation, 
or by lying within t he same cosmi c filament as 
RXJ1347 (see §4.2 in lLu et alJboiOl for a discus- 
sion). RXJ1347-SW lies outside both the Chandra 
and XMM-Newto n observations for t he cluster, so 
we use ROSAT (jVoges et al.lll999l ) to place an 
upper limit on its luminosity of Lx ^ 7 x lO^** 
erg s~^ w hich, using the Lx I Mggs.^m scaling re- 
lations in Zhang et al.l ( 201l[ ). predicts an upper 
limit to its mass (Mgas,500 ^ 10^"* M0) that is in 
agreem ent with the mass estimate from ILu et all 
(|2nin[ ). Ano ther nearby oyerden sity of galaxies 
was noted in Bradac et al.l ( 20081 ) immediately to 
the SW of RXJ1347, coincident with a mass exten- 
sion seen via weak-|-strong lensing. We also do not 
detect this object in the Chandra image and place 
an upper limit on its luminosity of Lx < 1.3 x 
lO'*^ erg s"^ at the distance of RXJ1347. 

There is some disagreement among the different 
lensing analyses as to the presence of the second 
massive dark matter halo that would be associ- 
at ed with eithe r cD2 or t he gaseous subcluster . 
In lBradac et~all (|2008l) and lMiranda et al.l (|2008h . 
their projected mass density maps show an ex- 
tension in the direction of the gaseous subcluster 
~30" (^2 70 kpc) southeast of the X-ray bright- 
ness peak. Miranda et al. (j2008l ) give an approxi- 
mate mass for this structure comparable with the 
primary cluster mass (^--^lO^^Mf;)) . On the other 
hand, the weak lensing analysis bv lLu et~aL ( 2010l ) 
shows no such structure, and does not even show 
an extension of the weak lensing signal in that di- 
rection. We discuss this further in terms of our 
proposed merger history for this cluster (see 

2.3. Radio Observations of RXJ1347 



Komatsu et al. I (ll999l) andlPointecouteau et al 



|l999, 2001) made the first measurements of the 
SZ effect in RXJ1347. The high an gular resolution 



(20" F WHM at 150 GHz) of the iKomatsu et al 



(|200l[ ) observations showed a strong enhancement 
of the SZ signal 20" southeast of the cluster cen- 
ter, at the same location as the surface bright- 
ness extension observed in X-rays (see Fig. [2] 
(a),(d)). Since the SZ intensity is proportional 
to the integrated gas pressure along the line of 
sight, this strong SZ signal was interpreted as be- 
ing due to hot gas, probably shock heated by the 



merger of the subcluster with the primary clus- 
ter. The presence of hot {kT > 15 keV) gas in the 
southeast region h a s been reported from Chandra 
(lAUen et all I2OO2I: iMiranda et al.l lioosh . XMM- 
Newt on dCitti et ahlbOo'Tb^ " and Suzaku observa- 
tions ( Ota et aLll2008 ). and is confirmed here. 
Recently. iMason et al. ( 2010[ ) and Korngut et al. 



()201l[ ) reported 90 GHz observations of RXJ1347 
with 10" angular resolution performed with the 
MUSTANG bolometer array on the Green Bank 
Telescope. In addition to confirming the enhanced 
SZ effect in the region of the subcluster merger, 
the MUSTANG observations (see Fig.HJd)) show 
a pronounced asymmetry in the projected pres- 
sure, with a high pressure ridge running north to 
south over ^ 0.5' (^197 h^^ kpc) approximately 
15" (~99 h-^ kpc) east of the X-ray peak. The 
southern end of this r idge marks th e peak of the 



SZ decrement, which Mason et al. ( 2010l ) inter 



pret as a shock between the primary cluster and 

infalhng subcluster (RX J1347-SE). 

We also note that iGitti et al.l (|2007al ) ob- 
served a radio mini-halo, consistent with the 
sloshing scenario (iMazzotta fc Giacintucci 2008t 
ZuHoneet al.ll201ll ). It is ^420 kpc in ex- 



tent and elongated to the southeast towards the 
location of RXJ1347-SE, with total flux den- 
sity 55.3±0.6 mjy centered on the cluster core. 
The highest resolution Ve ry Large Array im age 
(1.7" X 1.2" restoring beam: lGitti et al.ll2007al) also 
shows extended emission in the core of RXJ1347, 
with a diffuse (i.e. point-source-subtracted) radio 
fiux density of 25.2±0.3 mJy. 

3. Our New X-ray Data Analysis of RXJ1347: 
ACIS-I Data Reduction, Imaging and 
Photometry 

RXJ1347 was observed for a total of 80 ks over 
three Chandra observations (ObsIDs 506, 507, and 
3592). The first two observations were performed 
with ACIS-S and the third with ACIS-I. We re- 
processed the l e vel 1 event files as described in 



Vikhlinin et al 



(|2005l) Fl The background flare 
rejection was done using t he ratio of the 2.5-7 
keV and 9.0-12 keV bands (iHickox fc Markevitch 



20061 ) ■ excluding regions with point sources and 



^AU spectral extraction and spectral analysis were per- 
formed using CI AO y er. 4.3, the CALDB 4.4.1 and Sherpa 
l lFreeman et al.ll200ih . 



5 



cluster emission. This resulted in a combined 
~73 ks clean exposure, on which we performed 
the further image processing tasks described be- 
low. Among the products created in this initial 
reduction are background images and exposure 
maps which incorporate the telescope aspect solu- 
tion, vignetting, energy-dependent effective area, 
and position and energy-dependent detector ef- 
ficiency. The background-subtracted, exposure- 
corrected image used for all subsequent position 
and imaging analyses is shown in Fig. [1] Fig. [UJa) 
and (b) and Fig. HI We detected point sources 
using the wavelet detection algo rithm wvdecomp 



(|1998[ ). The 30 point- 



outlined in lVikhlinin et al. 
source regions within 2.5' of the X-ray peak were 
checked by eye and then excluded from all subse- 
quent analyses. 

Proper estimation of the X-ray background is 
critical for spectral analysis of the faint extended 
emission from clusters. To measure this, we used 
the blank sky data sets from the Chandra calibra- 
tion database (CALDB ver. 4.4.1). The aspect 
solution for our observation was applied to the 
blank sky exposures, remapping the background 
events to our source coordinate frame. We then 
performed an overall background normalization by 
measuring the ratio of the count rates at high en- 
ergies (9-12 keV, where the telescope's effective 
area is small) in our source images to that in our 
background images. For each of our spectral ex- 
traction regions, an identical region on the detec- 
tor was chosen from the blank sky image to pro- 
vide a background estimate. Associated ancillary 
response files and redistribution matrix files were 
created using CIAO ver. 4.3. 



3.1. Gas Properties Across the Surface 
Brightness Edges 

From the source, background and exposure 
maps, we create surface brightness images with 
different pixel binning: full resolution (no bin- 
ning) and block n with n x n pixel binning (where 
n = 2,4,8). The full resolution and block two 
images reveal two sharp edges east and west of 
the X-ray brightness peak (see Figs. [Hand [2][a)). 

In RXJ1347, the presence of two surface bright- 
ness edges to the east and west indicates that the 
core gas may be sloshing, however it is not com- 
pletely clear from which angle we are viewing it. 
A spiral pattern is characteristic of sloshing core 



gas however, depending on the degree to which 
the merger is inclined with respect to the plane of 
the sky, this spiral pattern may appear elongated 
along one axis. In RXJ1347, we may be viewing 
such an inclined spiral pattern elongated approxi- 
mately along the N-S direction. This can be seen 
in Fig. 1 beginning at a radius of ^ 10" away 
from the X-ray peak and position angle of 140° 
(measured counter clockwise relative to west), and 
turning counterclockwise to a radius of ~ 5" with 
position angle 52°. Such an elongation could in- 
dicate either a highly anisotropic ICM or gravita- 
tional potential into which the gas is expanding 
or, more likely, that the perturbing object came 
along a trajectory that was inclined with respect 
to the plane of the sky (see ii4.3l for a discussion). 

Rather than the spiral pattern, if what we are 
viewing is instead the concentric crescent-shaped 
edges associated with a merger primarily along our 
line of sight, we then face the difficulty in explain- 
ing how both cD2 and the gaseous subcluster could 
be undergoing a merger with the primary cluster, 
since their projected distances would then suggest 
real distances which extend far outside the cluster 
core. In this highly inclined merger scenario, we 
would also expect to see a significant difference in 
the redshifts of cD2 and cDl, which we do not.. 
We discuss the observational constraints and im- 
plications of both the line of sight and plane of the 
sky scenario in more detail in i i4.3l 

From the surface brightness images, radial sec- 
tors were chosen to encompass the east and west 
edges. We measure the radial surface brightness 
profiles within these sectors and fit them with a 
two-component model for the gas density with a 
power law inside the surface brightness edge and 
a /3-model outside with the form: 



ne(r) = ni 



for r < rj 
for r > r, 



(1) 

where rj is the radius of the surface brightness 
edge, 7i is the power law slope interior to the 
edge, Tc is the core radius of the /3-model, ni com- 
bines the model normalization with the emission 
measure normalization for each model in their re- 
spective regions, and J„^ is the magnitude of the 
discontinuous density jump at the edge. 

The square of this density model is integrated 
along the line of sight and compared to the 
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Fig. 2. — (a) Same image as in Fig.[T]with annular sectors overlaid for surface brightness and spectral extraction, (b) Same as 
in Fig. [T] with 2x2 pixel binning, showing the spatial extent of the annular regions along with the large scale emission from the 
cluster, (c) Temperature map produced as described in i|3.2l (d) MUSTANG 90 GHz SZ map (Korngut et al. 2011) showing 
the SZ decrement ridge along the eastern side of RXJ1347. Colors in (d) are of SZ decrement (lower value indicating stronger 
integrated pressure) measured lowest (black) to highest (white) with a logarithmic color stretch. We also highlight the shock 
front which shows as a hot region in the temperature map and an SZ depression (indicating high integrated pressure) in the 
SZ map. All images except (b) are shown at identical positions and scales. Color bar is linear scale in temperature going from 
cool (purple) to hot (white). A color version of this figure is available in the online journal. 



observed surface brightness profile using a chi- 
squared minimization. Due to the temperature 
dependence on the emissivity of the gas, and thus 
the surface brightness, this is an iterative process, 
where we first fit the surface brightness assuming 
a uniform temperature. Then, upon deprojecting 
the temperature profile using the resulting density 
model (see §3.3|) . we use the deprojected tempera- 
ture profile to again fit for the emissivity-corrected 
surface brightness. Convergence is obtained typi- 
cally after only one such iteration. 



We plot the resulting radial profile for the east- 
ern surface brightness edge (red points, labeled 
"east" in Fig. El^a)), and find that the surface 
brightness profile across this feature may be ad- 
equately represented by Eq. [T] with a jump of 
J„^=2.2±0.6 (see Table ^ at the radius of dis- 
continuity (r - 9.2" - 60 kpc). The den- 
sity model is plotted in Fig. [3t. As discussed in 
§1.11 given sufficient time since the initial slosh- 
ing disturbance, a surface brightness discontinuity 
on one side of the cluster should be accompanied 
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Fig. 3. — Radial quantities plotted for east and west regions shown in Fig. [2 a) and (b). (a) Radial surface brightness, (b) 
deprojected temperature, (c) electron density, and (d) pressure profiles for the regions east (solid x's, light (red) shading) and 
west (open dashed squares, dark (blue) shading). Spectral extraction and fitting are described in ^3.31 Shadings in (c) and (d) 
are from the 90% confidence intervals on the model parameters. A color version of this figure is available in the online journal. 



by a second disco ntinuity on the opposite side at 
a different radius ( Ascasibar fc Markevitcb 20061 
; and m.l\ in this work). Indeed, the western sur- 
face brightness jump hes ^ 4.6" (~ 30 kpc) 
from the X-ray peak, closer than the eastern edge 
which hes ~ 9.2" (~ 60 h'^ kpc) from the X- 
ray peak. Surface brightness and density plots are 
shown in Fig. ^a) and (c) (blue points, labeled 
"west"). That it is closer to the X-ray peak in- 
dicates that the western edge is the younger of 
the two, resulting from the most recent passage of 
the core gas across the cluster potential minimum. 
This feature is also well fit by a two component 
density model, with a power law inside the edge, 
a /3- model outside (Eq. ([1])) and jump magnitude 



2.0 



+0.6 
-0.5- 



3.2. Temperature Maps 

In Figure [2t we present a temperatu r e ma p 



for RXJ1347, also shown in lZuHone et al.l (|2010f ) 



Temperature maps reveal interesting structures 
and inhomogeneities in the cluster gas which 
might otherwise be overlooked by only examin- 
ing the X-ray surface brightness distribution. 
To generate the map , we use the method of 



Lap, 
1), 



Markevitcb et al. ( 2000D . in which a series of 6 



images in non-overlapping energy bands (0.5-1.5, 
1.5-2.3, 2.3-3.5, 3.5-4.5, 4.5-6.5, 6.5-9.5 keV) are 
created along with corresponding exposure maps. 
Each surface brightness image is then smoothed 
identically with a variable length scale Gaussian 
kernel, such that the smoothing scale at each pixel 
decreases roughly as the square root of the number 
of counts in that pixel (out to a brightness of ~0.02 
counts arcsec^, below which the smoothing scale 
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is constant). This method has the effect of apply- 
ing a smaller smoothing scale in areas of bright 
cluster emission and a larger smoothing scale in 
the low surface brightness regions, resulting in ap- 
proximately the same relative statistical accuracy 
across the interesting bright regions of the cluster. 

The surface brightness images are then used 
to fit a MEKAL model spectrum for each image 
pixel, binned to those same energy bands, and 
a chi-squared minimization is performed allowing 
only the model temperature and normalization to 
vary. Because this method does not assume any 
a priori shape to the spectral region, as opposed 
to tessellation or contour binning methods (e.g., 
[Sanders. 2006), we may be confident that any fea- 
tures shown, particularly near the center, reflect 
an accurate (projected) morphology of the cluster 
gas. 

The temperature map shows large-scale inho- 
mogeneities of the cluster gas, indicative of an un- 
relaxed cluster. Of particular note in this map are 
the cluster cool core, the presence of the eastern 
and western cold fronts seen in the surface bright- 
ness image, two very hot {kT >25 keV) regions, 
and the shock front between the primary and sub- 
clusters, also seen in the SZ map (Figure [2](d)). 

3.3. Spectral Extraction and Fitting 

After fitting the X-ray surface brightness across 
the edges, we extract spectra in these same sectors 
(see Fig.[2ja) and (b)). The spectral bins are nec- 
essarily wider than the surface brightness bins as 
we require at least 2000 source counts (background 
subtracted) in the 0.8-7.0 keV range per spec- 
tral bin to constrain a thermal model with small 
enough uncertainties to identify the temperature 
jumps predicted at the cold front edges. Spectral 
extraction and fitting were perfor med in an anal 



ogous way to our previous work IjJohnson et al 



20101) . using the CIAO suite of tools. Spectra 
in each spatial region were binned with a min- 
imum of 30 counts per energy channel and the 
models were fit over the range from 0.8 to 7.0 
keV. For each spectrum, we fit a thermal mekal 
model convolved with an intergalactic absorption 
model (xsmekal * xszwabs). The cluster's red- 
shift was fixed at the median redshift o f its clus- 
ter m embers (^=0.451) determined by ILu et al 



19901) at the cluster position. 

We account for the gas lying in projection along 
the line of sight by using a variation of the algo- 
rithms projct for XSpec and deproject for Sherpa. 
Our algorithm takes a series of radially extracted 
spectra, fits a model to the outermost bin first 
and then moves inward. For each inward step, an 
additional model is added to the fit, whose gas 
temperature is fixed at the values obtained from 
the outer bins. The difference between deproject, 
projct, and our method lies in our separate deriva- 
tion of the three-dimensional gas density rather 
than leaving the density as a free parameter in 
the model fit. The density for each radial bin is 
calculated as described in H3.1l and the ratio of the 
emission measure in the projected region to that 
in the emission region is used as a normalization 
constraint for each of the projected regions. We 
present these radial quantities in Table [T] 

It should be noted that due to the small ef- 
fective area of Chandra at high energies, it is 
difficult to constrain high temperatures (> 20 
keV) with confidence, even with high numbers 
of counts. This is evidenced in the 90% confi- 
dence limits on the temperatures in the hottest 
radial bins (see Table [1]). Despite this limita- 
tion, we do find that the eastern and western 
edges identified in the surface brightness images 
are indeed cold f ronts, as expected from t he ga s 
sloshing model ( Ascasibar fc MarkevitchI l2006f) . 
From the deprojected temperature and density 
profiles, we are able to calculate the pressure pro- 
file across each edge which, in the absence of high 
radial velocities, should be approxim ately contin- 
uous ( Ascasibar fc Markevitch l2006l) . In order to 
produce this profile we multiplied the mean tem- 
perature by the mean density in each spectral bin. 
The mean pressure profile does not appear very 
smooth, however this is due primarily to the large 
spectral bins required to constrain the tempera- 
tures, and the large uncertainties in those temper- 
atures. Deeper X-ray observations are necessary 
to decrease the uncertainties in the spatial temper- 
ature distribution. We note, however, that even 
with the coarse temperature resolution, including 
both the temperature and density uncertainties, 
each of the profiles is continuous at the 90% con- 
fidence level (see Fig. [3l[d)). 



(|2010l) with an absorbing hy drogen column den 
sity Nh = 4.85 x 10^" cm~^ (jPickev fc Lockman 
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Table 1 



Radial Temperatures for RXJ1347 



Region 


R 

["] 


kT 

[keV] 


X2 (dof) 


West 


70.0 
25.8 
9.5 
2.5 


15.4t|^ 

24.0+5'? 
6.4t?;? 


191.7 (197) 
230.4 (297) 
195.2 (257) 
87.0 (108) 


East 


95.2 
34.9 
13.7 
4.5 


15.0t^3;« 
26.8ti^«« 
12.51?'? 


94.0 (80) 
110.7 (136) 
191.2 (214) 
112.9 (160) 



Note. — MEKAL model temperatures for 
the regions shown in Fig. [2ja) and (b). 
(Col. 1) Spectral extraction region identifier. 
(Col. 2) emission- weighted radius of spectral 
bin. (Col. 3) deprojcctcd MEKAL model 
temperature. Errors on kT are the 90% confi- 
dence interval, Ax^=2.71, for one interesting 
parameter). (Col. 4) for model and num- 
ber of degrees of freedom. 



Table 2 

Edge Properties for RXJ1347 



Cluster ID 


Edge Radius ["] 


T2 
Ti 


"el 


Pi 
P-l 


East 
West 


4 6+°-2 


^■°°-0.40 


r, 99+0.60 
^•^^-0.58 


1 40+°-®'^ 
^■^"-0.50 

u-or_g 2g 



Note. — Values in each column are the ratios of various gas quantities 
(temperature, gas density and gas pressure) across the eastern and west- 
ern surface brightness edges in RXJ1347 (see Fig.[lJ, where A^i represents 
the value of the gas property A'^ inside the front (r < r^^g^) and N2 is the 
value outside the front (r > r^dge). The errors presented on the ratios 
are based on the 90% confidence limits from those quantities. 
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3.4. The Gaseous Subcluster 

The gaseous subcluster to the southeast of the 
X-ray core hes within the halo of RXJ1347 (see 
Fig. [T]), making a study of its spectral proper- 
ties problematic. We isolate its emission by se- 
lecting a circular region in the surface brightness 
image containing its emission (see Fig. 2]) along 
with three other identically sized regions at the 
same distance from, but at different position an- 
gles with respect to, the X-ray peak. The flux 
we measure in the region surrounding the sub- 
cluster includes some cluster emission, which we 
measure as fx,sub+ciust = 3.8 ± 0.5 x 10^^^ erg 
s~^ cm~^ and compare with the mean flux in the 
other three regions {Jx,other = 2.3 ± 0.4 x 10"^'^ 
erg s~^ cm~^)0 So our estimate of the subclus- 
ter flux is the difference of these, or fx. sub — 
1.5 ± 0.9 X 10~^^ erg s~^ cm~^. This subcluster 
flux translates to a subcluster luminosity (at the 
distance of RXJ1347) of Lx,sub = 1.1 ± 0.6 x lO''^ 
erg s~^. Usmg the L-Mgas scaling relations from 
Zhang et al.l pOllh l^l this corresponds to a sub- 
cluster gas mass of Mgas,sub ^ 4.8 ± 2.5 x 1O^^M0 
within r.r^ nn- Assuming a, core gas mass fraction 
/ga«=0.1 (jVikhlinin et al.ll200dh . this gives a mass 
for the subcluster of Msoo.sub = 4.8±2.5 x lO^M©. 
Given the merger scenario we propose in fjll the 
subcluster would have likely lost some portion of 
its mass during its first pericentric core passage. 
This suggests that our mass estimate is probably 
a lower limit. 



4. Discussion: The Merger History of 
RXJ1347 

While previous works agree that RXJ1347 ap- 
pears to have uiidergone a sigii ificant merger 



(|Gitti et al.l l2007at ILu et al.l mm . there is lit 



tie consensus on the details. The cluster has a 
complex morphology at nearly every wavelength, 
and there is no single merger scenario thus far 
presented which can explain the unique charac- 
teristics of this cluster. Specifically, the key ob- 
servables that must be explained by an acceptable 
merger model are as follows: 



^AU fluxes here computed in the 0.5-2.5 keV energy band. 

Errors are systematic variation of la based on the three 

regions surrounding the cluster. 
^ There are separate scaUng relations for both disturbed and 

undisturbed systems, and we use the former. 



• The presence of two sloshing edges: one at 
~ 9.2" (~60 kpc) east of the X-ray peak, 
and a second ~ 4.6" (~30 h^-^ kpc) to the 
west of the X-ray peak. 

• A very hot region located ^ 17" (~113 h^^ 
kpc) to the SE of the primary cluster. 

• A second cD galaxy with no detected X- 
ray gas atmosphere located ^ 18" (118 
kpc) directly east of the primary cluster and 
~ 11" (72 kpc) north of the gaseous 
subcluster. 

4.1. Previous Merger Scenario 

The presence of the high pressure ridge along 
the east side of the cluster (see Fig.[2](d)), connect- 
ing the second bright central galaxy (cD2 in Fig.[T]) 
to the gaseous subcluster, suggests t hat these ob- 
jects are all related. As noted in Im ason et al.l 
these features are both consistent with cD2 
having moved north (in projection) along the east 
side of RXJ1347, with the gaseous subcluster hav- 
ing been either stripped or shock heated (or both) 
as a result of this passage. This is the simplest 
model that agreed, until now, with all of the ob- 
servations. One problem with this model is that 
sloshing edges are produced well after the passage 
of the subcluster, which does not agree with the 
current location of either cD2 or the gaseous sub- 
cluster. 

4.2. Sloshing Constraints on the Merger 
History - A Second Crossing 

As discussed in simulations of gas slosh- 

ing predict a minimum time (~0.3 Gyr) should 
elapse between the first pericentric crossing {tper) 
of the perturbing cluster and the formation of the 
first sloshing cold front, irrespective of the merger 
axis inclination with respect to the line of sight. 
For the subclustei0 to have perturbed RXJ1347, 
it should now be a minimum distance away from 
pericenter for us to see sloshing cold fronts in the 
core. For example, if we assume a reasonable 
value for the subcluster's average velocity since 
tper, say 1000 km s~^, then it should be at least 



*Here we make special distinction between the gaseous sub- 
cluster, RXJ1347-SE, and the subcluster in general, com- 
prising cD2, the gaseous subcluster, and their combined 
dark matter halo. 
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Fig. 4. — Surface brightness image as in Fig.[Tl showing the regions used for estimating the subcluster emission. Each circular 
region is identical in radius, with their centers equidistant from the cluster X-ray peak (shown with vectors pointing from each 
circle). The source region used for the subcluster emission is labeled (blue, solid), whereas the other circles (blue-dashed) are 
used to represent the cluster emission at that same radius. A color version of this figure is available in the online journal. 



~300 kpc away in 0.3 Gyr. This minimum dis- 
tance is further constrained by the observation of 
a second sloshing cold front, which occurs only 
after the density peak has fallen past the poten- 
tial minimum, reached apocenter, and has then 
fallen back again0 Simulations show that this 
second front forms, at the earliest, ~0.6 Gyr af- 
ter iper, increasing the predicted separation dis- 
tance between pericenter and the current location 
of the subcluster to ~600 kpc. This is consider- 
ably larger than the observed (projected) separa- 
tion distance (^^160 kpc) between RXJ1347 and 
the location of cD2 (which should also be the lo- 
cation of the subcluster, as both the dark matter 
and the cD galaxy pass essentially collisionlessly 



This is assuming a near plane of the sky viewing angle. 



through the envelope of RXJ1347). This discrep- 
ancy indicates that the simple scenario where the 
subcluster is merging for the first time is not suffi- 
cient. We therefore propose a new scenario, which 
we will refer to as the "second crossing" scenario. 

4.3. A New Merger History For RXJ1347 

The simple interpretation of RXJ1347 now ex- 
periencing the first crossing of the subcluster does 
not correctly predict the presence of the slosh- 
ing cold fronts observed in the cluster core. It 
is also inconsistent with the simulation result that 
the primary cluster's core gas does not separate 
from its potential minimum until after the first 
passage of the subcluster. We propose instead 
that the subcluster first passed on the west side 
of RXJ1347, moving in a (projected) southerly di- 
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rection, initiated core gas sloshing upon reaching 
pericenter in its orbit, and has now returned for a 
second pass, this time from the south, and much 
closer to the core. During this second passage, 
the subcluster had its core gas stripped as it ap- 
proached pericenter, and is now near this point 
in its orbit. This "second crossing" scenario ad- 
dresses each of the three cluster observables noted 
in g 

The presence and location of two sloshing cold 
fronts around the primary cluster's core. By al- 
lowing the core gas sloshing to be initiated on 
the first subcluster passage, we solve the discrep- 
ancy between the minimum time since pericen- 
ter (~0.6 Gyr) and the current location of the 
subcluster. In this scenario, we require that the 
time between pericentric crossings of the subclus- 



diet is the case (cf. Tittlev & HenriksenI 


2005 


Ascasibar & Markevitchl2006; ZuHone et al. 


201C 


Roediger et al. 


.2011). Simulations also predict 



that the first cold front should be produced on the 
opposite side of the cluster from the perturber's 
pericentric location. For RXJ1347, we propose 
that the initial subcluster pericentric crossing was 
on the west side of the cluster. Thus the first 
cold front should be on the east side. Our mea- 
surements of the cold front radii (see Table [2]) 
indicate that the eastern cold front is farther from 
the X-ray peak (9.2") than is the western cold 
front (4.6"), i.e. it was formed earlier, in agree- 
ment with the "second crossing" scenario. We 
also note that the lowest mass subcluster that 
has been simu lated to produce slos hing is of or- 
der 10^3 Mq (jRoediger et al.ll201ll) . As shown 
in WSAl our estimate for the subcluster gas mass 
alone is Mgas ~ 10^^ Mq, putting its total mass 
Mt ~ 10^^ Mq. This implies the primary cluster 
is of order 10 times the mass of the subcluster, 
which is well above the mass ratio limi t required 
to initiate sloshing ( ZuHone et al. 2010l ). 

A very hot region located between the primary 
and sub- cluster, and the current location of cD2 
and the gaseous subcluster. Shocked gas between 
the primary cluster and the gaseous subcluster 
is predicted in all merger simulations for which 
the subcluster possesses gas of its own. This 
is because, being initially bound to the subclus- 
ter, the subcluster gas velocity exceeds the sound 
speed of the main cluster as it approaches peri- 



center in its orbit. The subcluster undergoing a 
second crossing will not change this. The shock 
from the subcluster's first crossing would have 
passed through the primary cluster's core already, 
likely contributing to the initial displacement of 
the R XJ1347's gas core fro m the potential mini- 



mum (jChurazov et al.l 120031 ). Before the subclus 



ter reached apocenter on its first pass, this shock 
would have traveled well beyond the cluster core, 
so that there would be little evidence of the front 
during the subcluster's second crossing. 

As the subcluster returned for its second cross- 
ing, this time from the south, its trajectory 
brought it closer to the core of RXJ1347, where 
it is now near pericenter, at ~ 18" (118 h^^ kpc) 
east of the primary cluster center. At this closer 
pericentric distance, the combination of the sub- 
cluster's high velocity and the ICM density of 
RXJ1347 increased the ram pressure on the sub- 
cluster's core gas enough to strip the gas from the 
dark matter halo. The gaseous subcluster is still 
traveling supersonically so we see a s hock front be- 
tween it and the core of RXJ1347. Bradac et al 



(j2008i see their Fig. 4) produce the highest reso- 
lution lensing map of the region around RXJ1347 
and show that the projected mass density is elon- 
gated towards the region in between cD2 and the 
gaseous subcluster. That the mass density elonga- 
tion does not lie directly on the gaseous subcluster, 
or on cD2, indicates that some fraction of the pro- 
jected mass density is lying north of the subcluster 
in the direction of cD2. This supports the notion 
cD2 has traveled along with the dark matter halo 
to the north of the gaseous subcluster. 

4.4. Simulated Merger 

A merger scenario similar to the one described 
above was i nodeled using the FLASH hydrody- 
namic code (jFrvxell et al.l I2OOOI) . The model ini- 



tial conditions are as described in IZuHond (l201lh 
for a different set of simulations with mass ratio 
i? = 10 and impact parameter b ~1 Mpc. When 
viewing the simulation in the XY plane (the plane 
of the clusters' merger axis), the subcluster ap- 
proaches on the western side of the main cluster 
on a southerly trajectory. Near pericenter, it ini- 
tiates sloshing in the primary cluster core. The 
subcluster reaches apocenter, then falls back to- 
wards the main cluster. Its proximity to the pri- 
mary cluster's core coupled with its relative veloc- 
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ity causes ram pressure on its leading edge that is 
strong enough to remove the subcluster gas core 
from the subcluster potential minimum. 

In Fig. [5] we show one particular snapshot in 
time of the simulated projected surface bright- 
ness, temperature and mass density of the cluster 
and subcluster and compare these with observa- 
tions of RXJ1347. Given the variety of distinct 
characteristics in RXJ1347 with which to com- 
pare, the qualitative similarities between the sim- 
ulations and observations are remarkable. In both 
the simulated and observed surface brightness im- 
ages (Fig.EJa) and (d)) the sloshing edges are ap- 
parent to the east and west of the primary clus- 
ter. The distances of the simulated fronts from 
the X-ray peak do appear different than those in 
the X-ray image. Specifically, in the simulated im- 
age, the eastern edge is closer to the X-ray peak, 
indicating that the sloshing there is at a more ad- 
vanced stage than in the observation. This is not 
surprising given that the current location of the 
density peak is dependent on many factors (peri- 
centric crossing time, impact parameter, mass ra- 
tio, the subcluster's relative velocity and orbital 
trajectory to name a few). The subcluster core 
also appears more concentrated in the simulation 
image compared to the Chandra image. 

Other merger related structures are visible 
when comparing the projected temperature maps 
(Fig. EJb) and (e)). Despite the poorer resolution 
of our temperature map relative to the simulation, 
we can still identify the location of the sloshing 
edges. We are also able to identify the cool region 
of the infalling subcluster. The peak in the SZ 
decrement is spatially coincident with the region 
just between the infalling subcluster and primary 
cluster, however the SZ signal extends further to 
the north, past cD2 (see Fig. EJd)). Our merger 
scenario accounts for the peak in the SZ decre- 
ment, since this corresponds to the merger front 
between the infalling subcluster gas and the pri- 
mary cluster's gas. At this front, we observe a 
temperature contrast and a pressure gradient (via 
the SZ signal) which are both consistent with it 
being a shock front. The significant SZ signal just 
to the north of where we identify this sho ck front 
(called the "ridge" in lMason et all (|2010[ )) is not 
explained by our merger scenario. It is possible 
that this ridge is associated with the eastern edge 
of the sloshing core of RXJ1347 and has no X-ray 



counterpart, however no such feature has been 
predicted in, or identified with, sloshing cores. 

In our temperature map (and i n other works; 
Ota et al. I I2OO8I: iMason et al.ll2010l) . there are two 
very hot (> 20 keV) regions which we do not see 
in the simulated temperature map. These hot re- 
gions may have been shock heated during the first 
passage of the subcluster. In the simulations, the 
entire cluster is heated as a result of this initial 
shock front, and by the second crossing the over- 
all cluster temperature has been raised from its 
pre-shocked state. The higher central tempera- 
ture in the simulated cluster core is evidence of 
this heating. The two hot regions, not associated 
with SZ decrements, which we find in our temper- 
ature map may indicate that the gas in these re- 
gions has cooled less relative to the simulation (i.e. 
less mixing), or that there has been some other 
mechanism (such as another shock front) which 
has maintained it at such a high temperature. 

Finally, we compare the simulation projected 
mass density to the locations of the two cD galax- 
ies (Fig. [51(c) and (f)). Because the galaxies and 
the dark matter in a merger should pass through 
nearly collisionlessly, the distribution of galaxies 
should trace the distribution of dark matter. In 
particular, since a cD galaxy resides at the bottom 
of its cluster's potential well, it should trace the 
location of the dark matter density peak. Com- 
paring the images, we see that the location of 
cD2 relative to cDl coincides with the location of 
the subcluster density peak relative to the main 
cluster density peak in the simulation. The loca- 
tion of the subcluster halo in the simulation also 
agrees with the elongation found in lens ing maps 



agrees witn tne elongation tound in iens ii 
(|Miranda et al.lboosi: iBradac et al.ll2008l) . 



5. Summary 

We have presented a new scenario for the 
merger history in the most luminous X-ray cluster, 
RXJ1347. This scenario is compared to optical. 
X-ray, and radio observations, and is able to re- 
produce several unique morphological features in 
the cluster. We perform a new analysis on the 
existing X-ray data that shows the core gas is 
sloshing, with sloshing cold fronts present to both 
the east and west of the cluster core. The mor- 
phology of the sloshing core gas places additional 
constraints on the merger history of the cluster: 
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Fig. 5. — Figure comparing a simulated cluster merger to observations of RXJ1347. On the left are the simulation snapshots 
with (a) projected surface brightness, (b) spectroscopic-like projected temperature and (c) projected mass density with linearly 
spaced contours of constant mass density. The observations on the right are (d) 0.5-2.0 keV surface brightness image smooth 
with Gaussian (FWHM=0.5"), (e) temperature map as in Fig. (2] (f) full resolution Chandra surface brightness image showing 
the location of cD galaxies. The spatial scales for each simulated image are the same, as are the spatial scales for each Chandra 
image. A color version of this figure is available in the online journal. 



namely that the interaction which initiated the gas 
sloshing occurred at least 0.6 Gyr ago, and that 
the perturbing object first passed along the east- 
ern side of the cluster. In order to reconcile the 
gas sloshing timescale with the current positions 



of the nearby gaseous subcluster and second cD 
galaxy, the perturbing subcluster should have first 
passed along the western side of the primary clus- 
ter, in a southerly direction and initiated sloshing 
near pericenter in its orbit. We are now seeing the 
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subclustcr when it is again near perieenter, this 
time heading north, where it has crossed close 
enough to the primary cluster's core that its gas 
has been ram pressure stripped. Wc make qual- 
itative comparisons between the multiwavelength 
observations for the cluster and high-resolution 
hydrodynamic simulations. Wc find that many of 
the unique merger structures around the cluster 
are well reproduced in the simulations. This par- 
ticularly well-studied cluster serves as an example 
of how core gas sloshing may be used to constrain 
the merger histories of galaxy clusters. 
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